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Abstract. Bayesian statistical models were developed for the 
number of tropical cyclones and the rate at which these cyclones 
became hurricanes in the North Atlantic, North and South Indian, 
and East and West Pacific Oceans. We find that there is small 
probability that the number of cyclones has increased in the past 
thirty years. The rate at which these storms become hurricanes ap- 
pears to be constant. The rate at which hurricanes evolve into cate- 
gory 4 and higher major storms does appear to have increased. We 
also investigate storm intensity by measuring the distribution of in- 
dividual storm lifetime in days, storm track length, and Emanuel's 
power dissiptation index. Wc find little evidence that, overall, the 
mean of the distribution of individual storm intensity is changing 
through time, but the variability of the distribution has increased. 
The cold tongue index and the North Atlantic oscillation index 
were found to be strongly associated with storm quality in the 
Western, and to a smaller extent, the Eastern Pacific oceans. The 
North Atlantic oscillation index was strongly associated with the 
increase in the rate of strong storms evolving. 



1. Introduction 



This paper carries on work we began in Briggd ( 2007I ). In that pa- 
per, we statistically characterized the number of tropical cyclones, the 
chance that these cyclones evolve into hurricanes, and the intensity of 
individual storms within a year. We asked whether or not the num- 
ber of storms were increasing, whether the chance that they evolved 
into hurricanes was increasing, and whether the year-by-year mean of 
intensity was increasing. We found that it was highly probable that 
the number of storms had increased, but that the other characteristics 
reamined constant through the 20th century for the North Atlantic. 

Here, we expand this work to include a global statistical model of the 
North Atlantic, North and South Indian, and West and East P acific 
oceans. Work on this topic is not new , e.g. IWebster et all (120051): an d 
e.g. for the Atlantic lEmanuell fl2005l ): IPielkd fl2005h : iLandseal fl2005h . 



Much interest naturally centers on wheth e r or not observed increases 
are due to global warminglTrenberthI (l2005[):lTrenberth and Sheal (120061 ) ; 
Landsea et al. ( 20061 ): Goldenberg et al. ( 2001 ). We do not seek to an- 
swer that question here. We only ask: have the number, and frequency 
of hurricanes and strong hurricanes that evolve from them, and the 
distribution of individual storm intensity changed since 1975? 

W e use the hurricane reanalysis dat abase (HURDAT) Jarvinen at 
al. (l2002h and, e.g. Landsea et al. (l2004 ). The remaining oceans 
data we got from the "best track" data on the UNISYS web site 
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(http: //weather. Unisys. com/hurricane/index.html ). This database con- 
tains six- hourly maximum sustained 1-minute winds at 10 m, central 
pressures, and position to the nearest 0.1° latitude and longitude from 
all known tropical storms from 1851-2006. A cyclone was classified as 
a "hurricane" if, at any time during its lifetime, the maximum wind- 
speed ever met or exceeded 65 knots. Obviously, this cutoff, though 
historical, is somewhat aribitrary and other numbers can be used: we 
also use the "category 4" or "super typhoon" classification of exceed- 
ing 114 knots. To investigate the realtionship of t ropical storms with 
ENSO, we also use the c old tongue index (CTI) iDeser and Wallace 
fll990h and e.g. ICail fl2003h. And we use the North Atlantic oscillation 
index (NAOI) from ljones et all fll997l ). 

There are obviously observational prob lems, even since 1975, with 
basins other than the Atlantic. FasuUo (120061 ) investigated observa- 
tional quality in the North Atlantic in the context of assessing trends 
and found that the data was reasonably consistent and that trends 
could be reliably estimated from it. Kossin et al. (120071 ) argue the 
opposite point and say that worldwide the data is rather inconsistent: 
they propose, and construct, a reanalysis of hurricane activity using 
satellite observations. Below, we will see that, particularly with the 
Indian oceans, the data is observationally biased through time. We do 
not seek to correct this bias, and remind our readers that all our results 
are conditional on the data we use as being correct. Where we suspect 
it is not, we give our best guess as to how to interpret the results. 

Section 2 lays out the statistical models and methods that we use. 
Section 3 contains the main results, and Section 4 presents some ideas 
for future research. 



2. Methods 



We once again adopt Bayesian statistical models. An important ad- 
vantage to these models is that we can make direct probability statet- 
ments about the results. We are also able to create more complicated 
and realistic models and solve them using the same numerical strategy; 
namely, Gibbs sampling. We do not go into depth about the partic- 
ular methods involved in forming or solving these models, as readers 
are likely familiar with these met hods nowadays . There are also many 
excellent references available, e.g. iGelman et al.l (120031 ). 



2.1. Number of storms. We suppose, in ocean i and year j of n, that 
the number of storms Sij is well approximated by a Poisson distribution 
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as m 

(1) ~ Poisson(Ai-,). 

We earlier found that this model well predicted the observed number 
of storms in t he Atlantic (s ee particularly Figs. 3-4 for a discussion 
of model fit in Briggsl . 2007 ). The number of storms in each ocean is 



assumed independent — this is backed up by the data where any corre- 
lation between the Sij for different i is low. We allow the possibility 
that Xij, the mean number of storms, changes linearly through time, 
and that th e CTI and NAOI may in f luence, or be associated with the 
mean CIT ( Eisner and Jagger . 2004 : Eisner et al. . 200 ll ). We use the 



generalized linear model 

(2) log(A,,) = /3^, + Pl.t, + /3^,fiTl, + /3|,NA0I, 

The superscipt s is to show we are in the number of storms portion of 
the model. The prior for each is 

(3) /3^j7fc,rfc~N(7fc,r,), A; = 0,1, 2, 3 

where is the precision (inverse of variance) , and to further abstract 
the mean and variance and thus allow them greater fiexability, we use 
the noninformative priors 

(4) ~ N(0, le - 6), ~ Gamma(0.001, 0.001). 

Now, it might be that the posterior of (31^ may be mostly or entirely 
positive for some but not all i; this would indicate that the mean 
number of storms is increasing in those oceans i. It may also be that 
the posterior of (31^, may be mostly or entirely negative for other oceans 
i'; this would indicate that the mean number of storms is decreasing in 
those oceans i'. All jS^^, over all oceans, are drawn from the distribution 
N(7f , Ti), and so if the mean of this distribution is greater than 0, then 
we can conclude that the mean number of storms is increasing across all 
oceans. Thus, we can look to the posterior of 7J to ascertain whether 
this is so. 

Once a tropical storm develops it, of course, has a chance to grow 
into a hurricane. If there are Sij tropical cyclones in ocean j and year i 
the number of hurricanes is constrained to be between and Sij. Thus, 
a reasonable model for the number of hurricanes hij in year j given Sij 
is 



(5) hij\sij, 9ij ~ Binomial (sjj, 9^ 



The parameter dij can be thought of as the proportion of hurricanes 
that develop from storms. It is possible, however, as with Ajj, that 6ij 
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is dependent on CTI and NAOI and that it changes through time. To 
investigate this, we adopt the following logistic regression model 

(6) log (y^) = + P'ut, + P^^CTl, + 

where the superscript h denotes the hurricane portion of the model 
and we again let the priors and hyperpriors be the same form as in the 
model for Sij. Like before, we will exmaine the posterior jSi^ across each 
ocean to see whether or not the it is likely that, for ocean i, the rate 
of hurricanes in increasing. And also if the posterior of 7^ has most of 
its probability above 0, we can conclude that the rate of hurricanes is 
increasing across all oceans. 

We carry the model one step further by asking about the chance 
that a category 4 or above (denoted by Cij) hurricane develops from 
an ordinary hurricane: the number of category 4+ storms is of course 
constrained to be between and hij. In analogy with the model above, 
we have 

(7) Cij\hij, ~ Binomial(/iij, ^tj). 

The parameter can be thought of as the proportion of category 4 
hurricanes that develop from ordinary hurricanes. And again, we allow 
^ij to be dependent on CTI, NAOI, and time. Thus, 

(8) log = Pl^ + f^uh + /92.CTI, + /?3^,NA0I, 

where the superscript c denotes the category 4+ portion of the model 
and we again let the priors and hyperpriors be the same form as in 
the model for Sij. Like before, we will examine the posterior j31^ across 
each ocean to see whether or not the it is likely that, for ocean the 
rate of major hurricanes is increasing. And if the posterior of has 
most of its probability above 0, we can conclude that the rate of major 
hurricanes is increasing across all oceans. 

The posteriors of each /^J*-, etc. are in logit or log odds space, so 
care must be taken in direct interpretation of any estimates. We leave 
them in this space so that a glace at the posterior density estimates 
tells the story: if most of the probabiltiy is above or below 0, we can be 
reasonably sure there is some effect of time. A simple way to transform 
these posteriors is by taking their exponentiation: the answers then are 
in odds and odds ratio space. We give some examples of this below. 

2.2. Measures of intensity. It may be that the frequency of storms 
and hurricanes remains unchanged through time, but that other char- 
acteristics of these storms have changed. One important characteristic 
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is intensity. We defi ne a three-dime n siona l measure of intensity, in fine 



with that defined in IWebster et aLl (120051 ): (1) the length m, in days 



that a storm lives; (2) the length of the track (km) of the storm over 
its lifetime; and (3) the power dissipation index as derived by Em- 
manuel, though here we apply this to each cyclone individually and do 
not derive a yearly summary. 

We approximate the number of days m to the nearest six-hours. 
Track length was estimated by computing the great circle distance 
between succesive six-hour observations of the cyclone, and summing 
these over the storm lifetime. The power dissipation index (PDI) is 
defined by 

I'm 

(9) PDI = / Vl^Jt 



where V^j^ax maximum sustained wind speed at 10m. Practically, 

we approximate the PDI — up to a constant — by summing the values 
(Kiax/100)3 at each six-hour observation. The PDI is a crude measure 
of the strength of the potential destructiveness of a tropical storm or 



hurricane, as cited by Emanuel (120051 ). Other than this measure, we 



say nothing directly about storm destructiveness (in terms of money 
etc.). 

It was found that log transforms of these variables made them much 
more managable in terms of statistical analysis. Transforming them led 
to all giving reasonable approximations of normal distributions; thus, 
standard methods are readily available. There is substantial correlation 
between these three measures, which we account of in the model below. 

First let, for ocean year j, and storm k (there are Sij storms in 
year j in ocean i), yijt = (log(m)ijfc, log(track length)^^^, log(PDI)ijfc)', 
i.e. a vector quantity. The index / will denote the Ith measure of y (i.e. 
Vijki = log(m)jjfc etc). Then we suppose that 

(10) log?/i,fc ~MVN(/ii,fc,A,,) 

i.e. a multivariate normal distribution where Aj^ is the 3x3 precision 
matrix for each ocean and year. We model the mean as before 

(11) fi^jki = P^u + Puitj + P^CTlj + Psu^AOlj, 1 = 1... 3 

where the superscript z denotes we are in the intensity portion of the 
model. We further let 

(12) N(7rH,</.w), r = 0,l,2,3 
and where these hyperparameters 

(13) 7rw~N(ari,6w) 
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where we use the noninformative priors ~ N(0, le— 6), bri ~ Gamma(0.001, 0.001) 
and 

(14) (pri ~ Gamma(0.001, 0.001). 

As in the models for number etc., we will exmaine the marginal pos- 
terior Pl'i across each ocean and dimension / = 1,2,3 to see whether 
or not the it is likely that, for ocean i, and dimension /, the mean in- 
tensity is increasing. Also analogous to the above models, since each 
(3^^i has mean vr^/, we can examine its posterior: if it has most of its 
probability above 0, we can conclude that the intensity is increasing 
across all oceans. 

Two priors were considered for the precision (inverse covariance) Aj^: 
a simple noninformative and a realistic, but more complicated one. The 
simple (and standard) prior assumes 

(15) Aij ~ Wishart(Jf,3). 

i.e. a Wishart distribution with t hree degress of freedo m and If is the 



3x3 identity matrix for ocean i fiGelman et al.l . l2003l ). 



Below, we give evidence that there is an unambiguous, probably 
linear, increase in the variance of intensity through time, where the 
increase is proportionally equal in each of its members. There is also 
substantial correlation between the parameters, which appears constant 
over time. The ratio of variances between members of intensity also 
appears constant in time. So we built an informative prior to take 
these observations into account. Let af indicate the variance of the Ith 
member of intensity, and let pim = p indicate the correlation, assumed 
constant across l,m = 1,2,3. Let dm = '^m/'^i- The variance was 
modelled as 

(16) 4i = co + ci(tj - 1990). 

where cq and ci represent the intercept and slope of the observed in- 
crease in variance (subtracting 1990 centers the time and helps elimi- 
nate computer roundoff error). Then afj2 = d20'fji and afj^ = d^afji- 
The covariance between elements 1 and 2 was modeled as p^/d^afji. 
the remainder of the covariances were handled in a similar manner. 
After the variance matrix Vij was formed, we let 

(17) Aij- -Wishart (V^ij, 3). 

In practice, with some notable exceptions mainly due to data in- 
tegrity, both priors for the precision gave about equal results, except 
that, as expected, the second model gave results (slightly) less certain 
than the first. The results were also not sensitive to exact values of 
p, Co etc. used. They are so insensitive, in fact, that the third prior 
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Vij = Ifj gave nearly identical results with ( fTTI) : the only difference 
between ( !T7|) and (ITSl) is that the former is allowed to change year by 
year (and over each ocean), and the later can only change over each 
ocean and is fixed in time. 



3. Results 

All co mputation s were carried out in the JAGS 0.97 Gibbs sampling 



software iPlummerl (120071 ) on a Fedora Core 6 platform. The first 2000 
simulations were considered "burn in" and were removed from the anal- 
ysis: 100,000 additional samples were calculated after this, with every 
10th simulation used to approximate the posterior distributions (the 
9 out of each 10 were discarded; this thinned the posterior simula- 
tions and helped to remove the small amount of autocorrelation of the 
simulations). 

3.1. Number of storms. Figure 1 shows the number of storms s from 
1975-2006 for each of the five ocean basins. There is no evident overall 
trend, though the number of storms seems to increase in some oceans, 
and decrease in others. Other plots and other formal statistics (not 
shown) also indicate that the number of storms, in any given year, are 
independent from ocean to ocean. 

The next picture shows h/s for each of the oceans. There is eviden- 
tally observational problems with the North and South Indian oceans: 
note the increase in the late 1970s of the ratio of hurricanes to storms. 
This is almost certainly due to improvements in observations and is 
not entirely due to natural causes. This being true, we nevertheless 
assume, for now, that the data is error free. 

The third picture shows c/ s for each of the oceans. Again, there is 
no evident overall trend, and the data quality of the Indian oceans is 
again suspect. 

The joint model for number of storms, and the ratio of evolved hur- 
ricanes and category 4+ storms was then run. Figs. |4}l6] and Tables 
[T}l3] summarize the results. 

Fig. m shows there is high probability that s has increased in the 
North Atlantic (NA) and West Pacific (WP), but decreased in the 
North Indian (NI), South Indian (SI), and East Pacific (EP). Overall, 
as registered by the posterior of the parameter 7, there has been no 
overall increase or decrease. This is backed up by the results from 
Table [H For the model of s, there is a greater than 97.5% chance that 
there was an increase in s in the NA, because all the percentiles, and 
in particular, the 2.5%-tile are greater than 0. The same is true of the 
WP, and the opposite is true for the NI and the SI. Overall, there is 
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Figure 1. The number of storms s five oceans from 
1975-2006. 



about equal probability on either side of 0, indicating that the best 
guess is to say "no change" in time. 

In f Briggsl . 200?! ) we went further and calculated the estimated prob- 
ability that the posterior of /5i was greater than for the NA. It was 
felt that the data quality of the NA was sufficiently good enough to 
state the results in such precise terms. We do not repeat these calcu- 
lations here, because doing so is likely to convey more certainty than 
is warranted. We repeat that the models we use are only as good as 
the data that go into them: the most precision we offer is in the form 
of the standard quantiles. 
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Figure 2. The rate h/ s at which tropical cyclones be- 
come hurricanes from 1975-2006. 



Fig. [5], and Table [U shows that there is some probability that the 
rate that hurricanes evolved from tropical storms has decreased in the 
NA, WP, and EP, though the result is not conculsive (there is less 
than a 97.5% chance that this is true). But there is good evidence 
(greater than 97.5% chance) that the rate has increased in the NI and 
SI. Whether this is due to natural causes or abberations in the data 
collecting mechanism, particularly in the SI, it is not possible to say. 
Overall, however, the evidence is nowhere near conclusive that the rate 
is increasing or decreasing worldwide. 
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Figure 3. The rate c/s at which tropical cyclones be- 
come category 4+ hurricanes from 1975-2006. 



The results for the model of c/h are more striking. It is still true 
that the observations for the SI are problematic, but, except for the 
NI, there is good evidence that the rate of category 4+ storms evolving 
from ordinary hurricanes has increased through time in the other oceans 
and overall. The odds of a category 4+ storm evolving increase, as 
estimated by the median overall column, by exp(0.04) = 1.04 times 
per year: over 10 years this is an increase in the odds of 1.5 times. 

Tables [2] and [3] show the influence that the CTI and the NAOI have on 
the different models. The results here are mixed, and ocean dependent, 
as might be expected. Increases in CTI are associated with a decrease 
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Figure 4. The posteriors of (3ii (the regressor for 
change in time) and 71 (the parameter for overall change) 
for each ocean i for the model of s. 



the number of storms in the NA, and to a slight extent in the NI, and 
it tends to be associated with an increase in the number of storms in 
the EP. Overall, there is not much effect. 

Similar results hold for the ratio of hurricanes to storms: increases in 
CTI tend to be associated with a decrease the rate at which hurricanes 
evolve in the NA and NI, and are associated with an increase with 
the rate in the WP. We will find that the CTI, as might be expected, 
plays a large role for storm quality in the WP. The SI and EP show 
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Table 1 . Common quantiles of the model parameter f3ii 
and 7i ( the parameters for change in time) for each ocean 
i. Recall that the interpretion for the models h/s and 
c/h are in terms of log odds space: these numbers should 
he exponentiated for odds ratios. Results ( increasing or 
decreasing) for which there is at least 95% certainty are 
highlighted in hold in this and all other Tahles. 
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no change; neither is there an overall effect. Increases in CTI tend to 
be associated with a decrease in the rate at which category 4+ storms 
evolve in the NA and SI, and are strongly associated with an increase 
in the rate in the WP and EP. Overall, the effect is small. 

The NAOI has little association on either the number of storms or 
the rate at which these evolve into hurricanes. But it does have a 
stronger association on rate at which category 4+ storms evolve: in 
every ocean, and overall, the effect is positive: an increase in a higher 
NAOI is associated with an increase in the rate at which strong hurri- 
canes evolve. 
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Figure 5. The posteriors of (the regressor for 
change in time) for each ocean i for the model of h/s. 



3.2. Measures of intensity. There was one missing track length (out 
of 1756) in the South Indian ocean. There were 28 (out of 622) missing 
PDIs in the North Indian and 106 (out of 1756) in the South Indian. 
The multivariate model we used requires that be no missing values. 
So we imputed the small number of missing values by forming linear 
regression models of log track (or log PDI) as a function of log m and 
log PDI (or log m and log track), and used these to make predictions 
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Figure 6. The posteriors of (the regressor for 
change in time) for each ocean i for the model of c/h. 



0.10 



at the missing values. We do not expect that this unduly effected any 
of the results. 

Figures [7H9] show the boxplot time series of individual storm (within 
a year and logged) m, track length, and PDI for each ocean. A simple 
regression trend line for the medians of these distributions is overlayed 
to help guide the eye as to the overall trend, if any. There is still the 
same data quality problem with the Indian oceans, but there are some 
trends evident in these pictures. 
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Table 2. The same as Table 1, hut for the regressor due 
to the CTI. 
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Within each year we estimated the covariance between intensity and 
plotted, in Fig. [TOl the estimated variance of log(m) (solid line), 
log(track) (dashed), and log(PDI) (dotted). Log(PDI) has been scaled 
by dividing by 4 so that all data fits on one picture. The is a clear in- 
crease in variance of each measure in Fig. [TOl similar across all oceans. 
Certainly, there is a lot of noise around this trend, particular (again) 
in the Indian oceans, but its existence is unambiguous. This represents 
one of our main findings. 

To illustrate this in a different way, we present Fig. [TTl which is four 
estimates of the year-by-year variance of intensity for the Western Pa- 
cific: the classical point estimate (solid line); posterior median based on 
the increasing variance prior (ITTI) (open circles); the posterior median 
based on the non-informative year-by-year prior; the posterior median 
based on the simple ocean prior (fT71) (dashed line). The posterior me- 
dians differ trivially for the priors (flTl) and the simplified year-by-year 
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Table 3. The same as Table 1, hut for the regressor due 
to the NAOI. 



Ocean 


2.5% 


50% 


97.5% 
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NA 


-0.16 


-0.04 


0.05 


NI 


-0.14 


-0.01 


0.11 


SI 


-0.13 


-0.05 


0.03 


WP 


-0.05 


0.02 


0.10 


EP 


-0.06 


0.03 


0.14 


Overall 


-0.11 


-0.01 


0.08 




h/s 






NA 


-0.31 


-0.08 


0.08 


NI 


-0.46 


-0.10 


0.07 


SI 


-0.14 


-0.01 


0.15 


WP 


-0.16 


-0.02 


0.12 


EP 


-0.24 


-0.06 


0.09 


Overall 


-0.24 


-0.06 


0.08 




c/h 






NA 


0.04 


0.26 


0.53 


NI 


-0.02 


0.25 


0.59 


SI 


0.06 


0.26 


0.47 


WP 


0.00 


0.21 


0.40 


WP 


0.03 


0.25 


0.48 


Overall 


0.05 


0.25 


0.46 



version. The prior (IT^ does not track very well with the other priors in 
this ocean, but it does better in other oceans, particularly the Indian 
(not shown). The Bayesian and classical estimates are in very close 
agreement for log(m) and log(track), but differ markedly for log(PDI): 
this is also found in the other oceans: the difference appears largest in 
this ocean, the ocean which also appeared to have the clearest increase 
in variance over time. This exhibits the typical Bayesian "shrinkage" 
that is found in multivariate estimation. Note: the classical (solid line) 
should certainly not be taken as the true value of the variance to which 
the Bayesian estimate aspires. It is merely another estimate, one that 
may be too high (just as the Bayesian estimate may be too low). In 
any case, all methods show a clear increase in variance through time, 
which is all that we are claiming holds. 

From the data in Fig. [TOl we estimated the parameters of the second 
precision prior (fTTl) . We estimate (by averaging the classical estimates 



CHANGES IN NUMBER AND INTENSITY OF WORLD-WIDE TROPICAL CYCLONES 

North Atlantic North Indian South Indian 




West Pacific 



East Pacific 



D) 
O 




1975 1985 1995 2005 



D) 
O 




1975 1985 1995 2005 



Figure 7. The time series of boxplots, for each year 
and each ocean, of log(m). A simple dashed trend hne 
of the medians in overlayed. 



over the oceans): p = 0.8 (.11), d2 = 2 (.56), = 6 (1.7), and 
Co = 0.3 (0.03), Ci = 0.004 (0.0026). The numbers in parentheses are 
the estimated standard deviations. The model results presented below 
are very insensitive to the exact values of these parameters. It turns out 
to be much more important whether the variance is allowed to change 
through time by each ocean, or that it is the same through time at 
each ocean. That is, we repeated all the analyses below, but this time 
with Vij = lij, i.e. nearly the same as the non-informative 'fiat' prior 
Vij = li, except we allow the variance to change year-by-year: the 
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Figure 8. The time series of boxplots, for each year 
and each ocean, of log(track length). A simple dashed 
trend line of the medians in overlayed. 



results below barely changed when this different, but eminently fair, 
prior was used. 

Figure fT2| and TablelH shows the posteriors of Pi^i and vr^i for log(m) 
for each variance prior. The posterior based on the increasing variance 
prior (fTTll is the solid line; the posterior based on the "flat" or constant 
non-informative prior (fT5l) is the dashed line. There is no difference 
in the NA: both priors indicate that there may be some evidence that 
log(m) has increased, but we cannot be very sure this is so. The same 
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Figure 9. The time series of boxplots, for each year 
and each ocean, of log(PDI). A simple dashed trend hne 
of the medians in overlayed. 



conclusion can be reached for the EP. No change is seen in the West 
and East Pacific. 

But for the North and South Indian, a complete opposite picture 
emerges (also partly seen in the West Pacific), particularly for the South 
Indian (where the data is most suspect). The posterior based on the 
flat prior (fTS!) gives good evidence that log(m) has decreased, but the 
posterior based on the increasing prior (fT7|) gives good evidence that 
log(m) has increased! This is partly because of the strong correlation 
between log(m) and log(PDI), and the strong rise over time of the later: 



CHANGES IN NUMBER AND INTENSITY OF WORLD-WIDE TROPICAL CYCLONES 

North Atlantic North Indian South Indian 




West Pacific 



East Pacific 



01 

o 
c 
ro 

CO 

> 



00 
o 



o 



o 



CM 
O 



O 
O 




T — I — I — I — I — I — r 

1975 1985 1995 2005 



O 

c 
to 

CO 

> 



00 
o 



CD 



CD 



O 




n — I — I — I — I — I — r 
1975 1985 1995 2005 



Figure 10. The estimated variance of log(m) (solid 
line), log(track) (dashed), and log(PDI) (dotted). 
Log(PDI) has been scaled by dividing by 4 so that all 
data fits on one picture. The is a clear increase in vari- 
ability of each measure. 



this is picked up in prior (|T7|) . We emphasize that this is not because of 
some bias in the informative increasing prior: a nearly identical picture 
arises if we instead use the year-by-year fiat prior. 

Now, the raw data picture of log(m) in the Indian oceans show de- 
creases, but there is a strong rise in the log(PDI) in the South Indian, 
which has about 3 times as many data points as does the North Indian 
ocean (and would thus receive more weight in the simulations). These 
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Figure 11. For the Western Pacific, four estimates of 
tlie year-by-year variance of intensity: tlie classical point 
estimate (solid line); posterior median based on the in- 
creasing variance prior f|T7j) (open circles); the posterior 
median based on the non-informative year-by-year prior; 
the posterior median based on the simple ocean prior 
(dashed hne). 



factors go a long way into explaining the difference in interpretation. 
In any case, any result for the Indian Oceans should be viewed with 
caution because these results are so sensitive on the prior used. 

The overall trend parameter for log(m) is near 0, meaning it is likley 
that log(m) has neither increased nor decreased across all oceans. 
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Figure 12. The marginal posteriors of /5i and tti for 
log(m) across all oceans. The solid line represents the 
posteriors with the increasing variance prior; the dashed 
line is for the simple, constant non-informative prior. As 
expected, in each ocean, the increasing variance prior 
produced wider — that is, less certain — posteriors. The 
results for both priors, however, are in rough agreement; 
but see the text for a discussion of the South Indian 
ocean. 



The same interpreation for log(track) and log(PDI) in Figs. [12] and 
[T3] as for log(m) can be made. The two variance priors give nearly 
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identical results, except for the Indian Oceans; and there is no overall 
trend. 

Table [5] shows how the CTI is related to intensity. For each dimen- 
sion of intensity, a higher CTI is associated with greater intensity in 
the WP, and to a smaller extent in the EP. To be clear: a high CTI 
means longer lived, longer travelling, and windier storms in the West- 
ern Pacific. Longer travelling storms in the NI are also associated with 
higher CTI. 
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Table 4. Common quantiles of the model parameter /Jn 
and 7i ( the regressor for change in time ) for each ocean 
i. Results are for both the increasing variance and non- 
informative ('flat') priors. There do not appear to he any 
overall trends. 
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50% 


97.5% 


2.5% 


50% 


97.5% 


log(m) 


NA 


-0.001 


0.005 


0.012 


-0.002 


0.004 


0.011 


NI 


0.000 


0.009 
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-0.014 


-0.005 


0.005 


SI 


0.004 


0.008 


0.013 


-0.016 


-0.011 


-0.006 


WP 


-0.002 


0.003 


0.008 • 


-0.009 


-0.005 


-0.001 


EP 


-0.004 


0.001 
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-0.007 


-0.002 
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-0.020 


0.005 


0.031 
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-0.003 
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-0.003 


0.006 


0.014 


NI 


0.000 


0.010 


0.024 


-0.006 
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0.012 


SI 
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-0.009 


-0.004 


0.001 


WP 


-0.003 
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-0.001 


0.004 


0.008 


EP 


-0.006 
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-0.004 
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Overall 


-0.021 


0.005 


0.031 


-0.024 


0.002 
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log(PDI) 








NA 


-0.004 


0.007 


0.019 


-0.006 


0.005 


0.016 


NI 


-0.002 


0.010 


0.022 


-0.008 


0.002 


0.012 


SI 


0.019 


0.027 


0.034 


0.007 


0.014 


0.021 


WP 


-0.006 


0.004 


0.014 


-0.007 


0.001 


0.009 


EP 


-0.011 


-0.001 


0.009 


-0.013 


-0.003 


0.006 


Overall 


-0.019 


0.009 


0.037 


-0.023 


0.004 


0.031 



The association of intensity with NAOI is presented in Table El 
There is little association in any ocean or overall, with the, perhaps 
surprising, exception that higher NAOI is associated with longer lived 
storms in the Western Pacific. 

4. Conclusions 

We find that there is good evidence that the number of tropical cy- 
clones over all the oceans basins considered here have neither increased 
nor decreased in the past thirty years: some oceans saw increases, oth- 
ers decreases or no changes. These results stands even after controlling 
for CTI and NAOI. These results are of course conditional on the model 
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Figure 14. The same as Fig. [12] but for log(PDI). 
Here, both variance priors give closer agreement except 
for the South Indian ocean. 



we used being adequate or at least it being a reasonable approximation 
to the data. As in our first paper, we make no predictions about future 
increases as it would be foolish to extrapolate the simple linear models 
we used into the future. 

We also found that the rate at which tropical cyclones become hurri- 
canes does not appear to be changing through time across oceans, nor 
is it much influenced by CTI or NAOI. The rate may be increasing in 
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Table 5. Common quantiles of the model parameter /3i2 
and 72 (the regressor for influence of CTI) for each ocean 
i. Results are for both the increasing variance and non- 
informative ('fiat') priors. 
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0.05 


SI 


-0.09 


-0.01 


0.07 
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0.06 


0.18 


-0.02 


0.09 


0.20 
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-0.10 


0.04 


0.18 


-0.07 


0.05 


0.17 



the Indian oceans, but it seems just as likely that flaws in the data 
would account for the results we have seen. 

There is good evidence that the mean rate at which major (category 4 
or above) storms evolve from ordinary hurricanes has increased through 
time, though the increase is small. 

We flnd almost no evidence that the mean of the distribution of 
individual storm intensity, measured by storm days, track length, or 
individual storm PDI, has changed (increased or decreased) through 
time over all the oceans. Again, there were certain noted increases in 
the Indian oceans, which may be real or may be due to flaws in the 
data: this is evidenced by the posteriors from these oceans being very 
sensitive to the priors used. We did, however, flnd an unambiguous 



eaiANGES IN NUMBER AND INTENSITY OF WORLD-WIDE TROPICAL CYCLONES 



Table 6. Common quantiles of the model parameter /3i3 
and 73 (the regressor for influence of NAOI) for each 
ocean i. Results are for both the increasing variance and 
non-informative ('flat') priors. 
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increase in the variance of the distribution of storm intensity over all 
oceans. 

It might be asked, if the overall number of storms has stayed con- 
stant, but the rate of category 4+ storms has increased, which implies 
that the number of category 4+ storms has increased, why has not, say, 
the mean of PDI (which is a direct function of wind speed) increased 
because of these stronger storms? This is most likely because the vari- 
ance of PDI has increased, meaning there has been an increase of both 
stronger and weaker storms, but that this change has balanced (in the 
sense that the mean has stayed the same). Of course, the variance of 
m and track length have also increased, leading to both shorter and 
longer lived, and shorter and longer travelling storms. 
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The CTI was particularly associated with storm quality in the West- 
ern Pacific: the rate at which hurricanes and category 4+ storms evolve 
increased with higher CTI. Higher NAOI also led, worldwide, to an in- 
crease in the rate at which category 4+ storms evolve. Higher CTI was 
also associated with longer and more powerful storms in the Western, 
and to a smaller extent, the Eastern Pacific oceans. Perhaps surprising, 
was that higher NAOI was associated with longer lived storms in the 
Western Pacific and nowhere else. 

In the Intr oduct ion we noted the hurricane reanalysis project of 
Kossin et al. (120071 ). The methods used here could certainly be used in 
this, and other datasets like it. These analysis from them would make 
a valuable comparison to the results presented here. 
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